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Abstract. Space provides unique opportunities to test gravitation. By using an interplanetary spacecraft 
as a test mass, it is possible to test General Relativity at the Solar System distance scale. This requires to 
compute accurately the trajectory of the spacecraft, a process which relies on radio tracking and is limited 
by the uncertainty on the spacecraft non-gravitational acceleration. 

The Gravity Advanced Package (GAP) is designed to measure the non-gravitational acceleration without 
bias. It is composed of an electrostatic accelerometer supplemented by a rotating stage. This article presents 
the instrument and its performances, and describes the method to make unbiased measurements. Finally, it 
addresses briefly the improvement brought by the instrument as far as orbit reconstruction is concerned. 
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1 Introduction 

With the ever-increasing precision of measurements, space has become a privileged place to test the two funda- 
mental theories which emerged during the 20 th century: General Relativity and Quantum Theory. In addition 
to providing a very clean environment, it opened new ways of testing these theories: as an example of interest for 
this article, precise navigation of in terplanetary spacecrafts allows probing the scale dependence of gravitation 

at the Solar System distance scale (jJaekel fc Revnaudll2006l) . 

Even if most experimental tests support General Relativity ( Will 2006() . there are still open windows for 



deviations. Indeed, the fact that these two fundamental theories are difficult to reconcile suggest that General 
Relativity may not be the final description of gravitation. The reason is that gravitation is the only interaction 
not having a quantum description. The validity of the Newton potentia l has been extensively tested fo r distances 
between the millimeter and the characteristic size of planetary orbits (Fischbach fc Talmadgelll999f). But there 



remain open windows outside this distance range for violations of the inverse square law (jAdelberger et al.ll2003l 
Fig. 4): below the millimeter or for distances of the order or larger than the Solar System characteristic size. 

Long range tests are performed using the motion of planets and interplanetary probes. Monitoring of 
the Moon and Mars delivers high precision tests o f the validity of General Relativity at these distances (e.g. 
Kolosnitsvn fe Melnikovf2004 : Williams et al.lll996h . However, t he navigation data of the Pioneer probes show a 



discrepancy with respect to the predictions of General Relativity ( Anderson et al.|[l998l 12002a : Lew et al"1 l2009). 



This discrepancy can be described as an anomalous acceleration directed toward the Sun with a roughly constant 
amplitude of approximate ly 8 x 10~ 10 m.s~ 2 . The origin of this anomaly is yet unexplained despite a huge effort of 
the scientific community ( Turvshev fc Toth 201fjl and references therein) : it may be an experimental artifact as 



well a s a hint of considerable importance for fundamental physics (jBrownstein fc Moffat 2006; Jack et fc Revnaud 



2005). At larger scales, the rotation curves of galaxies and the relation between redshift and luminosities of 



supernovae are accounted for by introducing respect ively "dark matter" and "dark energy", which represent 
25 % and 70 % of the energy content of the Universe ( Frieman et alj|2008l ). Since these dark components have 



been introduced on the basis of gravitational observations solely, the hypo thesis that General Relativity is not 
a correct description of gravitation at these scales needs to be considered ( Carroll et al. 20041) . 
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It is therefore essential to test gravitation at all distance scales. To this extend^ 



several mission con- 
cepts have been proposed t o improve the experiment made by the Pionee r probes (lAnderson et al. 2002b; 
Bertolami fc Paramosl[2007t Ijohann et aTll2008t IChristophe et al.ll2009l l201lt IWolf et al.ll2009t) . In many pro- 



posals, the addition of an accelerometer being able to measure without bias the non- gravitational acceleration 
of the spacecraft is central. ESA included this idea in the roadmap for fundamental physics in space (|ESAll2010h 
and recommended the development of accelerometer compatible with spacecraft t racking at the 10 pm .s~ 2 level. 



This article presents such an instrument, called the Gravity Advanced Package ( Lenoir et al. 2011bl ). First, a 



description of the i nstrument and its p erformances is given. Then, the method used to make absolute measure- 
ment is described ([Lenoir et al.ll2011aj) . Finally, the expected improvements of the orbit reconstruction process 
using the instrument are briefly discussed. 



2 The Gravity Advanced Package 

The Gravity Advanced Package is an important technological upgrade for future fundamental phys ics missions in 
space . It is composed of two subsyst ems: MicroSTAR is a three-axis electrostatic acc elerometer (jjosselin et al.l 
Il999h based on Onera's experience (jTouboul et al.lfl999l iTouboul fc Rodriguesll200lh . and the Bias Rejection 
System is a rotating stage with piezo-electric actuator used to rotate MicroSTAR around its x axis. The 
accelerometer aims at measuring the non-gravitational acceleration of the spacecraft but other quantities are also 
measured. In fact, MicroSTAR measures the components of the vector a on its three orthogonal measurement 
axis called x, y and z : 
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where ms and m A are the masses of the satellite and the proof mass respectively, F^^g is the non-gravitational 
force acting on the spacecraft, fl is the rotation vector of the instrument with respect to a Galilean reference 
frame, 1 is the vector between the center of mass of the satellite and the instrument (lever arm), FgL^ is the 
gravity of the spacecraft, and the last term in parenthesis is the gravity gradient expressed in term of acceleration, 
Fg t ^ s and F e G xi 



being the gravitati onal forces acting o n the satellite and the proof mass respectively. All 



these additional terms can be removed ( Lenoir et alj|201 lbh . 

Of course, the measurement is plagued by scale (<5fci„) and quadratic factors (&2i/), by bias (&„) and by noise 
(n„), such that the actual measurement on the axis v £ {x; y; z} is : 



= (1 + Ski u )a^ + k 2u a lJ 2 + b v + n v 



(2.2) 



where a v is the projection of a on the axis v. 
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Fig. 1. Left: Square-root of the power spectrum density of the accelerometer noise. Right: Square-root of the Allan 
variance of the accelerometer noise. The curves are plotted for different cut-off frequencies. The oscillations for 2nrf c < 1 
are not described by the simplified formula (|2.4|) . 



The measurement noise is characterized by the following power spectrum density (jLenoir et al 



2011bT) for a 
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measurement range of 1.8 x 10 4 m.s 2 (cf. fig. [!}: 
S n (f) = (5.7 x 1CT 11 m.s- 2 .Hz- 1/2 



4.2 mHz 



0.27 Hz 



(2.3) 



The characterization of the noise can also be given in term of Al lan varianc e A(t, f c ), where f c is the cut-off 
frequency. A simplified expression of the Allan variance (AVAR) (lAlfanlll966h for 27rr/ c 3> 1 is: 



An(T,f c ) = (5.7 x 10~ n m.s- 2 .Hz 



-1/2 



1 f 3 
— + 2 ln(2) x 4.2 mHz + -if-jr x (0.27 Hz)" 4 
2r 3tt 2 t 2 



(2.4) 



Figure [T] shows the dependence of Allan variance with respect to integration time and cut-off frequency without 
the 27tt/ c 3> 1 approximation. 

The Bias Rejection System rotates the accelerometer around the x axis of a monitored angle called 8. 
Assuming that the quadratic factors are equal to zero and that the scale factors are perfectly knowrQ, the 
quantities measured along the axis y and z are : 

m y = [cos(6>)ay + sm(6)az] + b y + n y (2.5a) 
m z = [— sin(0)ay + cos(8)az] + b z + n z (2.5b) 

with (jx G {Y; Z}) being the components of the acceleration in the reference frame of the spacecraft. 



3 Method to remove the bias of the electrostatic accelerometer from the measurements 

The method for removing the bias of the instrument consists in flipping MicroSTAR. The underlying idea is 
that when 8 = rad, the accelerometer measures the quantities m y = ay + b y and m z = az + b z , and when 
8 = 7r rad, it measures m y = —ay + b y and m z = —az + b z . Subtracting these measurements allows recovering 



T he complete method , which 
"" l2011al) . It is 



the external signal without bias, under the assumption that they are constants, 
can handle time variations of the external signal and of the bias, is described in ([Lenoir et al. 
shown in particular that the modulation signal, i.e. the time variation of 8, must fulfill some conditions in order 
to correctly remove the bias from the measurement (cf. eq. (|3.ip ). 

The modulation signal is supposed to be periodic, t being the period. Moreover, the measurements used 
for data post-processing are the ones made when the angle 8 is constant with time and only two positions are 
considered: rad and it rad (so that assuming k% v = is not restrictive). On the contrary, the measurements 
made when the accelerometer is rotating are not used because they may be spoiled by unwanted signal (vibration, 
fictitious acceleration). The duration of the rotation per period is called Tm and will be referred to as the masking 
duration. 

Assuming that the signal to measure and the bias are affine functions of time for each modulation period 
(which is correct if r is small compared to the characteristic variation time of the signal and the bias), the 
conditions for completely removing the bias from the measurements read 



r/2 f-r/2 i-r/2 

t m(t) cos(9(t))dt = t 2 m(t) cos(6(t))dt = 



r/2 



m(t) cos(8(t))dt 



r/2 



r/2 



(3.1) 



where m(t) is equal to when the accelerometer is rotating and 1 when it is not. These conditions allow deriving 
the time variation of the angle 8, which is shown in figure [2] (left). The pattern depends on the ratio of the 
masking duration Tm and the period r, as shown by figure [2] (right). 

With such a signal and after post-processing, it is possible to recover the mean of the external signal without 
bias over a modulation period r. It is possible to characterize these unbiased measurements in term of noise. 
The level of noise depends on the modulation period r and the masking time Tm but is approximately white 
as shown by figure [3] (left). For r = 600 s and Tm = 200 s, the uncertainty on the unbiased measurements is 
4.2 x 10 -12 m.s -2 (this value is the integral of the power spectrum density shown in fig. [3] (left)). This allows 
reaching a precision of 1 prn.s^ 2 for an integration time of three hours. 



* T hese assumptions are made for simplicity purpose. The complete treatment of the problem is presented in i Lenoir et al.1 
I2011af) . 
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Fig. 2. Left: Modulation signal (-) for a masking duration representing 33.33 % of the modulation period r = 1 arbitrary 
unit. Two periods are represented, separated by circles (o) Right: The figure describes how the modulation signal changes 
when the ratio Tm/t changes (one period is represented). The abscisse is time and the ordinate is Tm/t expressed in 
percentage. The curves indicate the values of the angle 9: between two curves of the same style 6 changes, between the 
curves (-) and (•) and between the curves (- -) and (- •) 9 = n rad, and 8 — rad elsewhere. 



4 Improvement of orbit reconstruction 



As mentioned in the introduction, the Gravity Advances Package (GAP) aims at improving the orbit reconstruc- 
tion of interplanetary probes so as to test the theories of gravitation. So far, models have been used to correct 
for the non-gravitational acceleration of the spacecrafts. But the computed orbit is then subject to errors due 
to uncertainties or inaccuracies in the models. By providing unbiased measurements of the non-gravitational 
acceleration in the orbit plane, the GAP enhances orbit reconstruction: it removes parameters to be fitted, it 
measures the fluctuation of the non-gravitational acceleration, and it removes correlations. 
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Fig. 3. Left: Square root of the power spectrum density of the unbiased measurements obtained after post-processing for 
a modulation period r = 600 s and a masking duration Tm = 200 s. Right: Square root of the power spectrum density 
of the solar radiation pressure noise (Fro hlich fc Le an 2004) in term of acceleration for a spacecraft with a ballistic 
coefficient of 0.1 m 2 .kg _1 at 10 AU of the Sun ()Biesbroekl 120081 ) . The bump at 3 mHz corresponds to the 5-minute 
oscillations of the Sun. 



There are several sources for the non-gravitational acceleration, the main ones being the direct solar radiation 
pressure and the anisotropic thermal radiation of the spacecraft. Figure [3] illustrates the interest of the GAP as 
far as the fluctuations of the non-gravitational acceleration are concerned. The power spectrum density of the 
unbiased measurement noise is compared to the noise on the radiation pressure expressed in term of acceleration. 
Whereas these fluctuations are not taken into account in the models, it shows that the GAP allows measuring 
them with a very high precision. 
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5 Conclusions 

The Gravity Advances Package, designed to improve the tests of gravitation at the Solar System distance scale, 
displays performances which will allow improving the accuracy of orbit reconstruction significantly. Indeed, with 
a carefully designed calibration signal for the rotating stage, it is possible to remove completely the bias of the 
electrostatic accelerometer and to obtain the mean non-gravitational acceleration of the spacecraft with no bias 
and a precision of 1 pm.s~ 2 for an integration time of three hours. This expected precision will allow reaching 
the 10 pm.s" 2 level for spacecraft tracking recommended by ESA. To do so, the accelerometer measurements 
will have to be taken into account during the orbit reconstruction process. 

The authors are grateful to CNES for its financial support. 
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